Potential Analysis manuscript No. 

(will be inserted by the editor) 



Harnack Inequality and Regularity for a Product of 
Symmetric Stable Process and Brownian Motion 



Deniz Karli 



Abstract In this paper, we consider a product of a symmetric stable process in R 
and a one-dimensional Brownian motion in R . Then we define a class of harmonic 
functions with respect to this product process. We show that bounded non-negative 
harmonic functions in the upper-half space satisfy Harnack inequality and prove 
that they are locally Holder continuous. We also argue a result on Littlewood- 
Paley functions which are obtained by the a-harmonic extension of an L p (R d ) 
function. 
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1 Introduction 

In the last 20 years, there has been an increasing interest in non-continuous 
stochastic processes. In particular symmetric a-stable processes play an impor- 
tant role in today's probability theory, and there has been a remarkable increase 
in the number of applications of symmetric a-stable processes. 
In this paper, we focus on a product process Xt which is the product of a d- 
dimensional symmetric a-stable process and a one-dimensional Brownian motion. 
We define a-harmonic functions with respect to this process in a probabilistic way 
and then study its applications. 
The organization of this paper is as follows: 

One of the most fundamental tools about harmonic functions is the Harnack in- 
equality. In section [3l we prove a Harnack inequality for a-harmonic functions in 
R d x R + . First, we show that the hitting probability of a Borel set with positive 
measure is positive. Then we prove the Harnack inequality using Krylov-Safonov's 
approach [6] . In section 2] we study regularity of these a-harmonic functions and 
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prove their Holder continuity. 

One of the ways to obtain an a-harmonic function is to extend a function on ~R d 
which can be taken as the boundary of K d x K + . This is also known as the Dirichlet 
problem for the upper-half space. We define the solution of this problem by means 
of the semigroup corresponding to Xt. Then in the last section, we state some 
results on Littlewood-Paley functions studied by P.A. Meyer and N.Varopoulos 
and then prove a result on a partial Littlewood-Paley function by using Harnack 
inequality. 



2 Preliminaries 

In this section, we give the details of the setup and introduce the notation which 
will be used for the rest of the paper. We consider the space R d with d-dimensional 
Lebesgue measure m and the upper-half space M. d x R + as a subspace of R d+1 . For 
simplicity, we identify the boundary of the upper- half plane, R d x {0}, with K d . 
We will denote the product measure m ® e a by m a for a 6 R + , where e a is point 
mass at a. Throughout the paper c, c±,C2,... will denote constants. Their values 
may vary from line to line. 

Now we introduce our process in details. Let Yt denote a d-dimensional right- 
continuous, symmetric a-stable process, that is, It is a right-continuous Markov 
process with independent and stationary increments and its characteristic func- 
tion is e - *'^ . It is known that Y± satisfies the scaling property with the factor 
of t 1 /™, that is, the processes (Y ct — Yq) and c 1 ^ a (Y t — Yq) have the same distri- 
bution for any c > 0. Let Pt denote the semigroup corresponding to Yt, that is, 
Pt(f)(x) = ~E x (f(Yt)). Here V x is a probability measure corresponding Yt started 
at the point x, and E x is the expectation with respect to V x . It is known that the 
measure P x (Y t 6 dy) is absolutely continuous with respect to Lebesgue measure. 
We will denote transition densities of this symmetric stable process by p(t,x,y). 
Unlike Brownian motion, there isn't any simple explicit formula for the transition 
densities. However, we will use the following estimate on p(t, x, y) : 

for some positive constants c\ and ci. Moreover, by [111 P.261] p(s,0,x) can also 
be expressed as 

roc 

/ (inu)- d / 2 e-W /(4ll) 9Q/2 ( S , W )d M (2) 
Jo 

where g a /2(s, u) is the density of an a/2 stable subordinator whose Laplace trans- 
form is given by J °° e~ Xv g a /2(s, v)dv = e~ sX 1 . Hence p(s, 0, x) is continuous and 
differentiable in the variable x. 

This stable process Yt forms the horizontal component of our product process 
if we think Mr as the horizontal and M + as the vertical component of x R + . 
On the vertical direction we let Z t be a one dimensional Brownian motion which 
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is independent from It, and so the product process is Xt = (Yt,Zt). We define 
the stopping time To = inf{£ > : Zt = 0}, which is the first time Xt hits the 
boundary which is R d x {0}. We note that To depends only on Z t by independence 
of Y t and Z t . 



In the classical context, the definition of a harmonic function can be given in 
both analytic and probabilistic way. We adopt the probabilistic definition and mod- 
ify it to define a-harmonic functions. In this paper, we call a continuous function 
u (x,t) on R d xR+ harmonic (a-harmonic) if the process u(Xt/\T ) is a martingale 
with respect to the filtration J- s = a{X r AT \T < s ) an d the probability measure 
pOM) f or an y starting point (x,t) G K d x R + . One natural way to obtain such a 
harmonic function is to start with a bounded Borel function / on R d and extend 
it to the upper-half space. One can write this extension by 

/>oo 

E (x '* ) /(^t )= / E x /(y s )P*(To G ds), (x,t) e R d x R+. (3) 
Jo 

Here, P*(To G ds) is the exit distribution of Brownian motion and its explicit form 
is well known, ft is the probability measure fit(ds) where 

Ht(ds) = J-e- t2 ^ s s-^ds (4) 

Zy/TT 

(see [7]). On the other hand, E x f(Y s ) can be written as P s (f)(x) = J f(y)p(s,x,y)dy. 
We note that this extension Q is a continuous function of (a;, t) (by (0]) and 
([2])). For simplicity, let's denote this extension by / as well and write ft(x) = 
f(x,t) = E^'V f(Y To ) for t > and f(x,0) = f(x). If we consider the process 
m( = f(X t /\T ) and the filtration T 3 = <j(X r/S T ', r < s) then it is easy to see that 
M( is a martingale under the probability measure p( x,t ) for any x and for any 
t > with respect to the given filtration. So the extension is a harmonic function. 

This harmonic extension can also be expressed as a convolution f(x,t) = f * 
qt(x) where qt{x) = J °° p(s,x,0)fJ,t(ds). If we define Q t by Q t = J °° P 3 (J,t(ds) 
then we can see that the extension of / is f(x,t) = Qt{f)(x). We note that Qt 
satisfies semi-group properties. 

In section [5l the main idea is based on the relation between a deterministic 
integral and a probabilistic integral. For the conversion between these two integral, 
we will use two tools, namely the Green function of Brownian motion and Levy 
system formula for the jump terms. If / is a positive Borel function, then the 
Green function for Brownian motion is given by 



E° 



To 



f(Z s )ds 



(z A a)f(z)dz 



for any a G R . The second tool, the Levy system formula, can be stated as follows. 



Theorem 1 Suppose f is a positive measurable function on R d xR d . Iff(x,y) 
zero on the diagonal then 



[ j f(Ys,Y s+ u)^ds 
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for any x G R . 

For the proof, we refer to the paper [3 . This property will be called as the Levy 
system formula throughout the rest of this paper. 

Now, let C be the infinitesimal generator of the symmetric a-stable process 
Y t on its domain T>(C), which is defined as lim t _,.Q+ (P t — I)/t. It can be given 
explicitly by 

= / [/(* + h) - f(x) - I { | h |<i } V/(i) • h]j^dh 

[TJ section 3.3]. We know that the semi-group Pt admits a square operator r which 
is given by the relation 

2r(f,g) = C{fg) - fC{g) - «?£(/). 

If we denote the point evaluation of r(f,g) at x by r x (f,g), then we have the 
following proposition. 

Proposition 1 If f is a bounded function inT>{C) and ft (x) = Qtf(x) as defined 
above, then 

r(f n-r f {Mx + h)-ft(x)} 2 
r.VU*)-cj |^p+^ dh. 

Proof According to the paper [71 Lemma 5 on p. 153], ft is in T>(£). Moreover, 
again by Theorem 2 on p. 147], ft is a l so in 

since the map x — > x 2 is in 
C 2 . By definition , 2_T(/ t ,/ t ) = jC(f 2 ) — 2f t C{ft). Now we need to consider two 
cases where 2 > a > 1 and 1 > a > separately. In the first case, we have 

ft(x)C(f t (x)) = f [ft(x)f t (x + h)- f?(x) - l m < iy ft(x)VMx) ■ h]rj-^dh 

and 

A/'M) = / ift{x + h) - f 2 (x) - ! { |^|< 1} V(/ t 2 )(x) • h)j^dh. 

So, 

2A(/ t , f t ) = f [ft(x + h)- f t ( x )\ 2 ^-dh 

+ I l m < 1} h.[2f t (x)Vf t (x)-V(f?)(x)}j^dh. 

Since V(/ t 2 ) = 2/tV/t, the second term is zero. Hence the result follows. In the 
second case where 1 > a > , we can drop the last term of the integrand of the 

generator, since the integral J TL{\h\<i} ^^ d+a dh is zero. Then a similar calculation 

yields to the desired result. 
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By Meyer P. 158], the quadratic variation of the martingale M/ is 

f tAT 

(M f )t = 2 / g{Y s ,Z s )ds (5) 
<9 



where g(x,t) = r x {f t Jt) + 



tAT„ 


2 



and r x (ft, ft) is as given above. This 



function appear in section [3] when we construct Littlewood-Paley functions. 

Before we end this section, let's list some short remarks which will be useful 
later. 

Remark 1 If m is Lebesgue measure, let P™ 1 " be the measure defined by 

l(x ' a) m(dx) 

and E'"" the integral taken with respect to the measure P m " . 
Remark 2 The law of Xt under the measure P ma is m. 

Remark 3 The semi-groups Pt and Qt are invariant under integration with respect 
to Lebesgue measure. If / is a bounded Borel function then 

J P t (f)(x)m(dx) = J f(x)m(dx) = j Q t (f)(x)m(dx). 
3 Harnack Inequality 

When one study harmonic functions on a domain, one of the most useful tools in 
harmonic analysis is the Harnack inequality. It allows us to compare values of a 
harmonic function inside a domain and it plays a crucial role for many applications. 
In this section, our main goal is to obtain a Harnack inequality (Theorem [3]) in 
this setup which was described in the previous section. To prove the inequality, we 
follow Krylov-Safonov approach (see [6]) and the method used by Bass-Levin in 
[3]. The idea is based on the fact that the hitting probability of a Borel set with 
positive measure is non-zero. So we focus on hitting probability first. 

Let's start by introducing the notation to be used in this section. Since we have 
different scaling factors in vertical and horizontal directions, we will consider the 
rectangular box D r (x,t) = PJ r (x) x B r (t), where 

D r (x) = {y e R :\xi — yi\< — — , i = 1, ...,d, x = (x lt ...,x d )} 

and 

B r (t) = {s e R : \s - t\ < ^}. 



For e G (0, 1), we define the box with e margin by Dr(x,t) = ~B%(x) x B( 1 _ e ) r (t) 
where 

-> £ d (i _ e 2/oL )r 2/a 

D r (x) = {y e R :\x t - yi\ < , i = 1, ...,d, x = (xi, ...,x d )}. 
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We denote the first hitting time and first exit time of a Borel set A by Ta and ta, 
respectively. If A is a point set then we will use T x instead of T{ x y. 

First we observe that expected exit time from the box D r (x,t) is comparable 
to r 2 . 

Lemma 1 Let e > 0. There exists a positive constant c such that 

E (y ' s W(*,t)) > cr 2 
for (y,s) G D e r {x,t) and D 2r (x,t) C K d x R+. 

Proof Let u > 0. Then using the independence of Y v and Z v first, and the scaling 
property second, we obtain 

E (j/,s W(*,t)) 

>^ (v ' s) {T Dr{Xit y,T Dr{x>t) > e 2 r 2 u) 
> e 2 r 2 uP^' s \r DAXtt) > e 2 r 2 u) 



> e 2 r 2 uP (^)( sup \ Yv - Y \ < eVa r Va ) sup | Z() _ Zq | < ^ 

= eVuP y ( sup |n-y l< ^ )P S ( sup \Z V -Z \<-) 

= eVuP°(sup \Y V \ < 1/2) P°(sup \Z V < 1/2). 

v<u v<u 

If we choose u small enough, then the last two probabilities above can be made 
bigger than 1/2. Then the result follows. 

Lemma 2 There exists a positive constant c such that 

E ( ^W(*,t)) < cr 2 

for any (y,s) G D r (x,t) and D 2r (x,t) C K d X R + . 

Proof By scaling we may consider the case r = 1, and so it is enough to show that 
E (2/ ' s) (Tr>i(z>*)) ^ c where D 2 (x,t) C R d x R+. Let S be the first time when Y v 
jumps of size larger than two. If we use the Levy system formula, 



SM <■ dh 



/•HAL r 

\s<i) = E ( y^ V i {lYv _ Yv _ l>2} =E ( y^ / / 

„<SA1 J ° J W> 2 



\h\ 



d+a 



dv 



dh 



h\>2 H d+C 



E M (SM) 



>cE ( ^ s) (S*Al; S> l) = cP {v ' 3 \S > l) = c[l-P ( ^ s) (S*< 1)]. 

Hence there is c' G (0, 1) such that p(«'")(S < 1) > c', and thus P^- S )(S > 1) < 
(1-c'). We note that T Dl{X}t) is smaller than 5. Thus P (y ' s) (T Dl(Xjt) > 1) < (1-c'). 
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Let's denote the usual shift operator by m . By the Markov property, 

P (y ' s W(*,t) > m + 1) = P (y ' s \T Dl{x!t) > m, T Dl(Xtt) o6 m >l) 

= E (v,«)(p(»,-)( TiMxji)O ro > l|JP m ); TDl(Xit) >m) 

= E (3/ ' s) (P x -"(r Cl(x>t) > l);r Cl(Xit) > m) 
< (l-c')P (y ' s) (r Dl(;Cit) >m). 

Then we obtain P (2/ ' s) (ri}i(<M) > m) < (1 - c') m by induction. This leads to 

OO 

E te ' s) (r Dl(x , t) ) < 1 + £ P ( ^ s) (r Dl(x!t) > m) = c < oo. 

m— 1 

Before we prove our theorem on hitting probabilities for general Borel sets, we 
consider the simple sets, which are the boxes of the form E x [a, b] C R d x R + 

Theorem 2 Suppose (y,s) G D2(x,t) and K = E X [a,b] is a rectangular box in 
Di(x,t) such that E C K d . Then 

P (y ' s) (T K < T D3(Xit) ) > c • m(£) • (6 - a) 

for some positive constant c and De(x,t) C R d x R + . 

Proof First, using the Levy system formula, we obtain 

\«<TD 3 (»,t) / 



To verify the last inequality, notice that \Y V — z\ is bounded above for any z G E 
under P^' s ). 
Next, we will show 



E (w.«) 



y 2{z„e[a,6]}rfwl > c(6- a). 



In this proof, we will use tJ^ and for the first exit time and the first hitting 
time of a Borel set A C R d by symmetric stable process Y t , respectively. Similarly, 

(2) (2) 

t a ; and ; denote the first exit time and the first hitting time of a Borel set 
A C R + by Brownian motion Z t , respectively. If A C R + is a point set then we 

(2) (2) 

will use the notation T^Jy = T„ ' . 

First, we note that T Da(Xtt) = A T ™ (t) = r-g^ A T^_ 3/2 t+3/2) . Set c a = 



( 3 2/« _2 2 / Q )/2. Then 



r Ca := inf{i; > : \Y V - Y \ > c a } < 



(x) 
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under p( y ' s ), since we have y £ z5 2 (a;). Hence 

e («m) ( f TD ^ 1{Zve[a b]} dv ^j > E ^ s > (^ r " AT3 1{z„e[a,6] } dv) 



f2) (2) 

where T3 = T t+3 /2 A ^t-3/ 2 - There are three possible cases: (i.) s < a, (ii.) s > b 
or (iii.) s £ [a, b]. First, assume that s < a and define a function h : R — »• R by 



if x < a, 

2 



(a - a) /2 if a; £ [a, 6), 

(6-a)(a;-^) if a; > b. 

Note that ft, > and h" = I( a ,b)i so that 

l E (o, s) ^ TC ° AT3 I {iJ „ 6[a , b]} = ^°> s \h{Z TcaA T 3 ) - h(Z )) 

> ^°' s \h{Z TcaA , 3 ) - h(Z );T t % /2 < T[% /2 A , 
= (h(t + 3/2) - h(s))V^ s \T^ 3/2 < T^% /2 A r c 

> (h(t + 3/2) - Hs^-^iT^ < T™ /2 A 
since s £ (t — 1, t + 1). Now we note that 

7i(t + 3/2) - fe(s) > (6 - a)(t + 3/2 - (b V s)) > (6 - a) 
since [a, 6] C (t— 1/2, t+1/2). Finally by translation invariance of Brownian motion 

P (0 ' t - 1) (^3/ 2 < T{% /2 A O = F (0 ' 1) (T^ < 2$ A O 

and so 

E ( °' s) (^ T °° ^ I{z„ 6 [a,6]} > c(6 - a). 
This proves the the first part. 

By symmetry, the case s > b follows from the same argument. If we take 
s' = It - s, a' = It - a and b' = It - b then 

E(0,s) (X°" Ar3 i{z ^ ia ' b]} dv ) = e(o,s,) {.C aAT3 i ^ b '^ dv ) 

by symmetry, and the last term is bounded below by c(d — b') = c(b — a). 

For the last case, assume s £ [a, 6]. We use the same function h(x) as in the 
case s < a. Arguing as in the case s < a, 

EiV ' S) {fo^ 1(z ^ um dv ^ ~ E(0,S) ^ f^^ 3 ljz ^ um dv 



{Z v £[a,b}} aV ) - ^ yj ^{Z v e[a,b]} 

2E {0 - 3 \h(Z TcaAf3 )-h(Z )) 

t+3/2 ^ " i t-3/2 



> 2 (ft(t + 3/2) - h(s))F^\T™ < l£> A r c 
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Moreover, by translation invariance of Brownian motion 

^ S) (T t % /2 < 4% /2 Ar c J > P (0 ' s) (Tj 2) 3 < Ti!\Ar c J = p(°- 2 )(T 5 (2) < if V C J. 
Since s G [a, b], we have h(s) < (b — a) 2 /2 and 

h(t + 3/2) - > (6 - a)(t + 3/2 - ^) - (b ~ a) 

> (&-a)(t + 3/2 -6) 

> (6 -a). 

Hence, 



e (im) 
So in any case 

E (v,«) 



(^ rD3( °' t) l{z„ e[ a, b] } <fo) > P (0 ' 2 >(T 5 (2) < T[ 2) Ar c J(6-a). 

(y l{Z„G[a,b]} <*« I > C(6- a). 



So we proved the hitting probability of a rectangular box is positive. We can 
extend this to any compact set in Di(x,t) using Krylov-Safonov's method which 
is based on covering the compact set with rectangular boxes. From here on, we 
denote the Lebesgue measure of a set A in R d x K+ by 

Corollary 1 There exists a non- decreasing function ip : (0, 1) — > (0, 1) such that 
if A is a compact set inside D\(x,t) such that \A\ > and (y,s) G D2(x,t) then 

P (y ' s) (T A <r A3(x , t) )>^(|^|) 

where D 6 (x,t) C R d X R+. 

Proof First define 

ip(e) = inf{P (z ' u) (T B < Tn 3 ( Z0 , U0 )) ■ (*0,«o) € R d x R + , («,«) G Da(*o.«o), (6) 
> e |Z>i(ao, wo)|, S C Di(zo, wo), B is compact, Dq(zo, uq) C ]R d x R + }. 
Set 50 = inf £• We claim that go = 0. Suppose not, that is, let's assume that go 

V (e)>0 

is strictly positive. This will give us a contradiction. First we note that go < 1. So 
we can find g such that (g + g 2 )/2 < go < g. Let 77 = (g — q )/2 so that g — ij = 
(<7 + <7 2 )/2- Let p > which will be chosen later. There is D 6 (z ,u ) C K d x K+, 
(z,u) G D2(zo, uq), and -B G l?i(zo,Mo) such that 

\B\ 

1 > in \T > 9 - »? (7) 

and 

P (z ' ll) (T B <r I , 3(zo , Uo) )<p.^(g) 2 . (8) 

Without loss of generality we may drop (zo,uo) and denote Di(zo,uo) by Di for 
i = 1,2,3. Now we construct the rectangular region D which is described in the 

\B\ 

Proposition (2) 7.2]. Since Bcfli and q > j — -7 = \B\, there exists a rectangular 

|-Di| 

region D which satisfies the following conditions: 
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1. D = UiRi where Ri is a cube and Ri denotes the cube with the same center 
as Ri but side lengths three times as long, 

2. the interiors of RiS are pairwise disjoint, 

3. \B\ < q- |Dn£>i|, and 

4. \BC\Ri\ > q - \Ri\ for all i. 

Since \B\ > q — r\ and \B\ < q ■ \D n Di\, we have 

lDnDll >W > !^l = I+l >q . 

q q 2 

Set D = D n D\, This results in \D\ > q. By the definition of p , we obtain 

^ U \T £) <T D3 )>p{q). (9) 
To obtain a contradiction we want to show that if (V, it') € D then 

vL*>«'\T B <T Da )>p-<p{q). (10) 
If this is true then by using the strong Markov property 
p(*.«)(r B < TDa ) > pC*.«)( Tfi < T B < r Da ) 

> E (z ' u) (P* T J5(T B < r C3 ); < r Da ) 
>p.^(g).p( z ^(T 6 <r I , 3 ) 



> P • ¥>(?) 



2 



This contradicts (JH1> - To obtain ([10]), first set F = ftnfli. Note that F is a 
rectangular box in D±. So Fi C -D3. Moreover Fi l~l B = Ri n B, since B G -Di. So 
we have |F, D B\ > q ■ \F t \, and if (z , u) £ F z then 



i( z '.«') 



(T Fl nB < Tp.) > 



by the definition of tp. 

If 0', u') G £> then (Y, u') eftflDi for some We can find a cube K cR d 
and A" x [c, d] C F. By the Theorem [2] there is p such that 

P (z ' X) (T Fl < r Ds ) > P (z> ' 3 (T Kx[Cid] < t D3 ) > p. 
Finally, the strong Markov property implies that if (z ,u') G D, then 

P («',«')( Tfl < TDa) > P (-'.«')(T Fi < r D3 , 7> <nB < 73) 
?(*'>«') f in>0'>«')/ 



= Hv* > (P^ z J (T Fl n B < r Fj | T F J; T Fi < rD 3 
= E^'' u ') (p Xr n (T FinB < r Fi );T Fi < rD 3 ) 
= ^( 9 )-P (z ' V) (T Fi <r J D 3 ) 

= p-<p(q)- 

This completes the proof. 
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Lemma 3 Suppose H is a function which is bounded, non-negative on R X K + 
and supported in (D2r{x,t)) c where Di r {x,t) C R d x K + . // (y,s), {y',s') e 
Dr(x,t) then 

E^H(X TuAXit) ) < c^y'^H(X TDr(Xit) ). 

Proof It suffices to prove the statement for H = lp with F C (D2 r (x, t)) c . Since 
Z v is continuous , H(X TD ) is non-zero if and only if Y jumps from D r (x,t) into 

(D2r(x, t)) c . Since D r (x,t) = ~3 r (x) x B r (t), it is enough to consider the sets 
F = E x [c,d] with [c, d] C B r (t) and the function H of the form J? = 1 F . Then 

»<"AT Dl , (ljt ) 

Since £ Z^ r (a;) and h £ E, \h — Y v \ is comparable to |/i — a;|, that is, there are 
constants ci and C2 so that 



ci 



/• dh f dh f 

J E \h-x\*+° J e \h-Y v \*+" - C2 J E 



dh 



E \h-x\ d + a ' 



So 



^(VS) f UATD r^,t) f dh f dh ^( VS ), 



and similarly 

^<vs) [ uAtd ^ r dh r dh ^ (vs)/ 



Now if we let u —> oo, and use the fact that cr 2 < E^' 5 - 1 (rc r ( Xit )) < cV 2 , then 

Observe that the above inequality also holds under P^ y ' s K Hence the result fol- 
lows. 

So far we have showed that any compact set with a positive measure is visited 
by Xt with positive probability given that the starting point is close. To prove 
the Harnack inequality, we will use this fact by defining a compact set on which 
the harmonic function takes large values. Since this set is visited with positive 
probability, one can define a sequence on which the function is unbounded. 

To simplify our notation, we fix a point temporarily and we will denote the 
usual rectangular box around this center by D r . 
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Theorem 3 There exists c > such that if h is non-negative and bounded on 
R d X R + , harmonic in Die and D 32 C K d X R + then 

h{y,t) < ch(y',t'), (y,t), {y',t') e D x . 

Proof By taking a constant multiple of h, we may assume that inf^ h = 1/2. Then 
there is (yo,to) £ D\ such that h(yo, to) < 1. We will show that h is bounded above 
in D± by a constant not depending on h. We show that if h takes a very large value 
in D±, then we can find a sequence in D 2 for which h is unbounded. 
Let (x, s) £ D\ such that h(x, s) = K. K will be chosen later and we can take 
K large. Let e > be small. By the Lemma Q, there is C2 such that if (x,s) £ 
l d xR + , r > 0, D r denotes the rectangular box D r (x, s) about the center (x, s), D% 
denotes the rectangular box with e margin D^.(x,s), H is bounded, non-negative 
and supported in D 2r then for any (y,t), (y ,t') £ D% 

E (y',t') H (x TDr ) < C2 E ( - y ^H(X TDr ). (11) 

By Theorem ©, 

F (y°^ { T D2r/3<TB J>c 3 \D 2r/3 \. (12) 
By Corollary ([T]), there is C4 such that if C C D r / 3 is a compact set with 
' > i then 



I A73I " 3' 

P (y ' t} (T c <T Dr )> Ci (y,t)eD 2r/3 . (13) 

Now let 77 = — , and £ = — A (c^r/). Let r be so that |-D 2r /3| = ^ . Note 

3 3 c 3 C4 QK 

that in this case r~ +1 = C5 K -1 , and so r = C5 ^ _a /( 2d + a ). Here we note that 

by taking a large enough value for K, we keep r small (that is r < ^) and so the 

rectangular boxes stay in the upper half space. 

Let A' = {(u,v) £ D r/3 : h{u,v) > C,K}. We claim that 

\A'\ K l 



\D r/3 \ ~ 2 

If not, then there is a compact subset A such that 

\A\ 1 



-D r /3 1 3 

Then (O and (13]) imply that 



1> h{yo,to) >E {v "' to) [h(X 




^ Ta< t d 1b 


} 


> (KE (yn > to) 




13 [Ta < r Bm 




> (KE {yo ' to) 




13 [T A < r Dr ] 


\ T D 2r/3 < r 5i J 


> (Kc 4 ¥ {yo - tn) \r D2r 


/a < T 5 J 





(14) 



> (Kc 3 a \ D 2r / 3 
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This contradiction shows that 



1A731 - 2 



Let C be a compact set in D r / 3 — A' such that 

_J£L>I. 

\D r/3 \ ~ 3 

Let H = l Dl h. If E {x ' s \h(X TDr );X TDr <£ D 2r ] > f]K then for any (y,t) G D r/3 

h(y,t) = E^[h(X TDr )] > E^[h(X TDr );X TDr # D 2r ] = E^[H(X TDr )] 
> c^E^iHiX^)} > c 2 \K > CK. 

But this is a contradiction, since A' is not all of D r / 3 . So we must have 

E {x ' 3) [h(X TDr );X TDr # D 2r ] < vK. 

Let's denote the supremum of ft on D 2r by M. Then 

K = h(x,s) =E {x ^[h(X TcATDr )} 

= E {x ' s) [h{X Tc );T c < T Dr ]+E^ s) [h{X TDr ); T c > r Dr ,X TDr G D 2r ] 
+ E {x ' s) [h(X TDr ); T c > T Dr ,X TDr $ D 2r ] 

< CKP (x ' s) [T c < r Dr ] +MP (x ' s) [T Dr < T c ] +vK 

< C KP (x ' s) [T c < r Dr ] + M (1 - P (x ' s) [T c < T Dr ]) + V K. 

So 

M > 1 - 7? - CP ( "' a) [Tc <T flr ] 
AT ~ 1 -P(z,*)[T c < rcj 

Then there is a ft > such that M > (1 + 2/3) iC, and, as a result, there is 
(x, s') G D 2r such that h(x' , s') > (1 + 0) K. 

Now suppose there is {x\, s±) G D± with h(xi, s±) = K\. We can then find n 
as above. And there is (#2, S2) G £>2n (#1, si) with ft(a;2,S2) = K 2 > (1 + f3) Ki. 
By induction we can create a sequence {{x%, Sj)} and corresponding {Ki}, {ri} so 
that {xi+\, Sj+i) G Z^r; (xi, Si) and i£j > (1 + j3) 1 K\. Note that 

^2\x i+ i - x t \ < c e K± 2/< - 2d+a) and ^ |sj+i - Sj\ < 07 K~ a / < - 2d+a \ 

i i 

So HKi > (2c 6 ) d+ f and Ki > (2c 7 )~ +1 then (a; i ,s i )'s are in D 2 , and ft(a;i, Si) > 
(l+/3) I_1 Ki which contradicts the fact that ft- is bounded. Hence Ki < (2c6) d+ ~V 
(2c 7 )^ +1 = c, and 

sup ft < c. 

61 
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4 Regularity Results 

One of the important application of the Harnack inequality is the regularity for 
the solution of elliptic PDEs. Bass and Levin [3] develop some techniques to prove 
Holder continuity and regularity results of some integral operators. They applied 
these techniques in the case of a jump Markov process whose kernel is comparable 
to that of a symmetric stable process. We will follow their method and show 
regularity results for our A-resolvent. 
We call the operator 

P OO P oo 

U x f(x,t)=E {x ' t) e- Xs f(X s )ds= e- Xs Q s f(x,t)ds 
Jo Jo 

as the A-resolvent of /. Resolvent plays an important role in applications. We 
study some properties of this resolvent. 

Theorem 4 If f is a bounded function on R d xM. + and it is harmonic on Di(x,t) C 
R d x K + ; then f is Holder continuous in D\(x,t), that is, there exists 7 > such 
that 

\f(v, s) ~ f(y', s')\ < c||/||oo|(y, a) - (y', s')\\ (y, a), (y\ a') G Di(x, t). (15) 

Proof Without loss of generality by taking a constant multiple, we may assume 
that ||/||oo = 1. First, we pick two points (y, s) , (y , s') G D\(x,t). Observing that 
D2(y,s) C Di(x,t), f is harmonic on D2(y,s) by the hypothesis. It is enough to 
consider the case when (y' , s') G Di(y, s). Because, otherwise, \(y, s) — {y' , s')\ > 1 
and hence 

\f(y,s)-f(y',s')\<2<2\(y,s)-(y',sT- 
By corollary [T] there is Co > such that for any compact set A C D r / 3 (y, s) with 
\A\ > \D r/3 (y,s)\/3 and r < 1/3 (so that D 2r (y,s) C M d x R+) 

p(w.«)(r A <T Dr{y<s) ) >co. (16) 

We fix co and choose ft G (0,1) close enough to 1 so that /3 2 > (1 — Co/4) > 0. 
Arguing as in lemma [3j if we take H = J[K£i X R+-r> r(y,s)] such that r > 2r then 
for u > 

P {v ' s) (Xr Dr(v .^u $ D r ,(y,s)) = E^ [H(X TDrMAu )] 

By Levy system formula and the fact that \Y V — h\ is bounded below for h $ 
D r '(y,s) and Y v G Z^ r (y,s), the last expression is bounded by 

^ / |T * +a dh dv < c(r') " 2 E (y ' s) [r DrM A «] . 

Let u — y 00. By lemma [2] we obtain 

P (tf, *H*T flp( ,,. ) A»*JV (»,*))<<* (^) 2 . (17) 
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We note that this inequality is valid even if r' is large and D r >(y,s) <£_ R d x R + , 
since the event considers only the jumps of Xt from D r (y, s) to K d x R + — D r ' (y, s) 
and a jump may occur only in the horizontal direction. Now let 

and denote infimum and supremum over nested rectangular boxes by 
<2fc = inf / and b k = sup /. 

D e k(y,s) D gk (y,s) 

We will show that b k — a k < f3 k by induction. First of all, it is clear that the base 
step bo — ao < /3 holds. Assume this inequality also holds for any i < k, that is 
bi — a% < f3 l . Now let A' = {/ < (dfc+i + 6fc + i)/2} n D 9 k+i (y, s). We may assume 
> \Dgk+i(y, s)\/2. Otherwise we can work with 1 — / instead of /. Take a 
compact set A C A' so that |A| > \Dgk+i(y,s)\/3. 

Let e > and pick (z,u), (z ,u) G Dgk+i(y, s) so that 

f(z,u) < a k+1 + e and f(z',u')>b k+1 -e. 
Since / is harmonic in D2(y, s), 

f(z,u) - /(*',«') = E<*>"> [f(X TDek{y:e) ) - /(*',«')] 

= /l +/2+I3- 

where 

= E (z ' m) [/(X ta ) -/(«', u');ta <r 0efc( „, s) ] 
h = E (z ' u) [/(X TD(t( „ f) )- /(/,«'); ta > r Dgkiv , s) ,X TDeklytS) e £> e »-i(y,a)] 

oo 

»=1 

X ^ efc («-) e D ek -i-i(y,s) - D g k-i(y,s)^ . 
We will bound each term in terms of powers of /3. First we note 

t s- ( ak + 1 + ^fc+l \ m>(«.«) ^ "\ 

/i < ^ 2 afc+1 J P l TA < rz V(^) J 

^ 2 J V A TD .*(»'-) J 

<^ P (-) (r A <r I , #fc („,.,) 

and 

J 2 < (6 fc _i - a fc _i) P (z ' u) (ta > r^ fc(y , s) ) 
<^ 1 P M (^>^ fc( „ s) ). 
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Hence 

/i+/ 2 </? fc - 1 



\^ {Z ' U) (ta < T Df)kM ) + 1 - P (2 ' u) (r A < r DthM ) 



<P k -\l-^)- (18) 

For the third integral, we will use ([TTJ) , the fact that 9 2 < f3/2 so that (1 — 
(^ 2 /7)) _1 < 2 and 6 2 < (/3c )/(8ci). Then 

OO 

i = l 

OO 

< J2 P k ~ l ~* ^ {z ' u) (X TD fc(H , 3) G D B k-i-i(y, s) - D gk -,(y, s) 



2 

nfc-1- 



i=l 



k-2 n 2 



< 2ci/3 

< (co/4)^" 1 . 



So by CLS) and (fT9|). 

f(z,u)-f(z\u')<p k -\l- ^)<P k+1 . 
This shows that b k+1 - a k+1 < f3 k+1 . 

Finally, fix k so that (y',s') G D e k(y,s) — D gk +i{y,s). We observe that 

\(y,s) - (y',s')\ > (8 k+1 A [8 k+1 ] 2 / a ) = [8 k+1 ] 2 / a 
since (y',s') g" D g k+i(y, s) and 2/a > 1. Hence 

log \(y, s) - (y', s')\ >(k + l)(2/a) log(0). 

Then 

< e (a/2) log(/3) log \(y,s)-(y',s')\/ log(0) 

where 7 = (a/2) log(/3)/ log(#). This is the desired result. 

Having the Holder continuity of harmonic function, we can discuss the conti- 
nuity of the resolvent. As we defined before, A-resolvent of / is given by 



U x f(x,t)=E {x ' t) e- Xs f(X s )ds. 
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We will show that 

\U x f(y,s) - U x f(y',s')\ < c||/|U(|(l/,«) - (y',a')\ A 1)^. 
First, we recall the resolvent equation 

(j3 — X)U\Up = U\ — Up. 



(20) 



(21) 



It is easy to verify this equation due to the semi-group property of Q s . This 
resolvent inequality reduces (|20[) to the case A = 0. So we will prove the following 
theorem first, and then we prove (|20[) by using this theorem. 



Theorem 5 Suppose f is a bounded function on R d x R + with compact support. 
Then there exists 76 (0, 1) such that 

\U f(y,s) - U f(y',s')\ < c(||/||oo + \\Uof\\ao) (\(y,s) - (y',s')\ A l) 7 

Proof First of all, without loss of generality we may assume D4(y, s) C R d x R + 
by scaling property. We observe that if \{y,s) — (y',s')\ > 1, the result follows 
immediately. So we assume \(y,s) — (y',s')\ < 1, that is (y',s') £ D±(y 7 s). By 
definition 

U f(y,s)=E^ f(X t )dt 



which can be written as 



U f(y,s) =E 



(y. s ) 



/ /(X t )dt + E^ / f{X. 



)dt. 



r < 1 will be chosen later. (We only assume r < 1.) By strong Markov property, 
the second term is equal to 



E 



Similarly, 



f(X t+TDAys) )dt = E 



(y.s)]g(yis) 



f(Xt+r DrM )dt 



(v, s ) 



Uof{y',s')=^y'^ 



E 



[Uof(X TDrM )). 



f(X t )dt + E^'^ [U f(X TDr(y ^)] 



Hence 



\U f(y,s) - U f(y',s')\ < ||/||oo [E {v - 3) r Dr{y ^ + E^' s \ DrM \ + 

\e^ [Uof(X TDr{y:B) )]-E^'^ [u f(x TDrltltt) )]\. 

Note that the first term is bounded by cr 2 1 1/| |oo- We also note that the function 
(z, u) — > E( z ' u > [Uof(X TDr(y > )] is harmonic in D r (y, s). So it is Holder continuous 
by previous theorem. Then the second term on the right is bounded by 



ll^0/||c 



\(y,s) - (y',s')\ 

r A r 2 / a 
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after scaling. Now let r = \(y,s) — (y' , s')\ a ^ 4 . Since \{y,s) — (y',s')\ < 1, we have 
\(y,s) - (y', S ')| 1/2 < \(v,s) - (y',s')| Q/4 , and hence 

\U f(y,s) - U f(y',s')\ < c \\f\U(y, s) - ( y ' , s')\ a ' 2 + c' \ \U f\ U (y, s) - (y',s')\^ 2 

< c" (H/lloo + ||C/o/||oo)|(y, S ) - (y', S ')| (Q ^ )/2 . 
Theorem 6 Suppose / is bounded and A > 0. Then we have 

\U x f(y,s) - U x f(y\s')\ < c|l/||oo(|(y, S ) - (y',s')\ A 1)^ 

Proof It is enough to show this for positive and compactly supported functions. 
Suppose / > and / has compact support. Define the function g by 



g{y, s) = f(y, s) - \U\f{y, s). 



(22) 



Note that 



XUoU x f = U f - U x f 
by resolvent equation (|21[) and hence 

Uog = U f - XUoUxf = Uxf. 
These equations can be verified by direct calculation. Using this equality, we obtain 

||E/0fll|oo = \\Uxf\\oo < c\\f\\oo 

and 

IMloo < I I /| |oo + A||[/a/||oo < C 1 1/| |oo. 

Finally, by previous theorem 

\U x f(y,s) - Uxf(y',s')\ = \U g{y,s) - U g(y',s')\ 

< c (IMloo + ||L%||oo) (1 A | (y, s) - (y\s')\y 
<c\\f\\ oa (lA\(y,s)-(y' ) s')\y. 



5 Littlewood-Paley Functions 

This short section is about Littlewood-Paley functions obtained from the a-harmonic 
extension of a function / on l d . After defining G-functions in this context we state 
some earlier results by P. A. Meyer |7:, 8 f 9, and N. Varopoulos [12_ and then we prove 
a result on a partial G-function which is close to the area functional used in the 
classical context. Let's denote the general Littlewood-Paley function by Gj and 
define it by 



Gfix) 



r 

Jo 



tg(x,t)dt 



1/2 



[f t ( X + h) - f t (x)f 
\ h \d+ a 



dhdt + 



f 

Jo 



dt 



dt 



1/2 
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where ft{x) = Qtf(x) as defined in the section of preliminaries. As one can easily 
see, the integrand g(x,t) is the square function which was defined before in ([5|. It 
has two components: one corresponding to the vertical process (Brownian motion) 
and the other one corresponding to the horizontal process (symmetric a-stable 
process). Let's denote these components by and G/, respectively. Explicitly, 
these two functions are 



[ftjx + h) - f t (x)}' 

\h\ d + a 



and 



GUx) 



f 

./() 



(It 



dh dt 



1/2 



1/2 



Some results on L p -norms of these G-functions are known. P.A. Meyer worked 
in [8] with symmetric Markov processes and proved an L p inequality for the case 



p > 2. When applied to our setup, we can state it as ||G/| 



< c 



for p > 2. 



On the other hand N. Varopoulos showed in his work [12] that this inequality 
can be extended to p > 1 for the Brownian component, that is |G^|| P < c||/|| p 
for p > 1. However the extension of this inequality for the general G-function 
is not possible which is pointed by M. Silverstein [10] . Here we discuss the part 
of the horizontal component on a parabolic-like domain. Although the extension 
of this L p inequality is not true for ~G* f when p G (1,2), we can obtain a partial 
result considering the operator only inside a parabolic-like domain. In the classical 
context, the area functional is defined in a similar way on a cone with a vertex 
at a point x G M. d . Since we have different scaling factors on each component, we 
consider a modification of this domain and study the part of G-function on the 
set {{x + h,t) £R d x I 



\h\ < t 2 / a }. For this purpose, define ~G* 



G f , a (x) 



(f t (x + h) - f t {x)Y 



dh 



\h\ d +L 



■ dt 



1/2 



First we note that the harmonic extension can be expressed as a convolution of 
the function with an approximate identity. To see this, we observe that transi- 
tion densities of the symmetric stable process Y t satisfy the relation p(st 2 , 0, y) = 
t~ 2d / a p(s,0,yt~ 2 ^ a ) by the scaling property. Hence we can write 



ft{x) 



f(x - y) <j> t (y) dy = f* <j> t {x) 



where 



and 



<j>(x) 



p(s,0,x)/j,i(ds) 



M^) = t~ 2d/a ^(x/t 2/a ) 



So ft{x) < cA4(f)(x) [H section 2.1] where -M(/) is the Hardy-Littlewood max- 
imal function. To see this, it is enough to note that (f> is radially decreasing and 
its Z^-norm equals one. Since the transition density p(s, 0, x) is obtained from the 
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characteristic function e s ' x ' by the inverse Fourier transform, it follows trivially 
that p(s,0,x) is a radial symmetric function. Moreover, by [11, P.261] p(s,0,x) 
can also be expressed as J °° (47ru) _<i/ ' 2 e — ' x ' g a / 2 (s,u)du where g a /2{s,u) is 
the density of an a/2 stable subordinator whose Laplace transform is given by 

J °° e~ Xv g a / 2 (s,v)dv = e _sA . Hence p(s,0,x) is radially decreasing in the vari- 
able x. The fact that <fi is radially decreasing follows from the previous line. 
It is also known that 



\\M(f)\\ p < c\\f\\ p , when feL p and p G (l,oo). 
Lemma 4 // / G L p , then for some £ = £(x, h, t) between ft(x) and ft(x + h), 

wm > c f jf°°t f e- 2 (f(x+h) - Mx)) 2 jj^dtdx. 

Proof Suppose / is positive and consider the map F(x) = (x + e) p for some 
e > 0. Clearly, F G C 2 (R + ). So we can apply the Ito formula using the martingale 
M t = f(X tA T ), and we obtain 

t 



(M/ + ef = (Af* + e) p +p (M/_ + tf^dMl 

Jo 



+ 



+ Yl [( M * + e ) p - ( M /- + e ) p - p( M *- + e ) p_1 ( M / - M L 



By [3 P.168], there is a positive function j(s) such that d((M^) c ) s > j(s)ds, and 
hence the third term on the right hand side is positive. Moreover, the convexity of 
the function F implies that the jump terms are also positive. So taking expectations 
first, and letting t — > oo and e — > 0, we have 

e^((m4) p ) • >:"• ( y. [wir - (m/_) p - P {M f s _y-\M! - m/_)] ) . 

\ s < T o / 

We recall that \\f\\ p = E m "((M^ o ) p ) [Remarks (HJ and ©]. If we denote by 
A(f, x, y, z) the expression 

A(f, x, y, z) = f p (x, z) - f p (y, z) - pf p -\y, z)(f(x, z) - f(y, z)) 

then we can write 

11/11? >E m » ( J2 [A(f,Y.,Y a -,Z a )]\ . 

Next, we will use the Levy system formula, invariance of P t under the Lebesgue 
measure m, and Green's function, respectively, and the right hand side of the last 
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inequality becomes 
E" 1 



" (/ Q T ° / A (f,Y s +h,Y s ,Z s )j^ds 



f E a (£° J A(f, y + h,y, Z.)j^ ds^j dy 



{a At) I A(f,y + h,y,t) TT ^j^dtdy. 



| ft | 



If we let a — > oo, then 



f\\ P P > t A(f,y + h,y,t) 



dh 

h\ d + a 



dt dy. 



If we use the Taylor expansion of x — > x p , then for some £ between ft{x) and 
ft (x + ft) we have 

A{f, y + h,y,t) = P^^le~ 2 (ft(y + h)- f t (y)f, 



and using this equality we obtain the result 

r r oo p 

\\f\\l>c t e~ 2 (f(x+h)-f t (x)y 



\h\ d + a 



dt dx. 



Theorem 7 If p e (1,2) and f e L p (R d ) then 

0/,a||p<c||/||p. 

Proof Without loss of generality we may assume that / > e > 0. Then it can 
be generalized to / € L p . By definition of the harmonic extension, we also have 
ft > e. By Lemma there is £ = £,(t,h,x) between ft{x) and ft(x + ft) such 
that 



\\f\\i>c t e~ 2 (f(x+h)-Mx)y 



dh 



\h\ d + a 



dt dx. 



Using this £, we can write 

Gf, a (x) = 



f (f t{x + h )- Mx )f^- ( ;, 

J{\h\<t 2 /°>} \ n \ 



d+a 



1/2 



ft/ 

Jo J{\h 



\<t'/"} 



(f t ( x + h)-f t (x)) 2 jj^dt 



1/2 



Temporarily, we fix t > 0. Let Rt be the rectangular box centered at (x,t) with 
side-length (^) 2/q , (^) 2/q , (^) 2/ °, ^, that is R t = D t/32 (x,t). Then R 32t C 
M d x R + . Let /3(a;, ft, i) denotes the linear path from (x, i) to (z + ft, i) where |ft| < 
t 2/a . This path can be covered by ra-many horizontal translations of the box R t , 
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say R\,R 2 t , R", such that R{ n R{ +1 £ for j e {1, 2, , n - 1} and i? t n i# ^ 0. 
Note that n can be chosen so that it does not depend on t, x, or h. Choose points 
from each pairwise intersection, say (x%,t) G Rt fl i?t, (a;2,t) 6 Jit PI i?t, ■•• , 
(x n) i) € Rt _1 H By using the Harnack inequality, we obtain 



/t(a; + /i) < cft(x n ) < c ft(x n -i) < ... < c ft(xi) < c ft{x). 

Hence, if |ft| < t 2/a , then f t (x + h) < c n+1 f t (x), and § < [/ t (x) V f t (x + h)) < 
c n+1 f t (x). This implies that £ 2-p < cf?~ p (x) < c' [M{f){x)] 2 - p . So G f , a {x) is 
bounded above by 



c' [M(f)(x)\^/ 2 



JO J {\h\<t 2 / a } 



e- 2 {f t {x + h)-f t {x)fj^<n 



\h\ d +t 



1/2 



By using the Holder inequality with 2/(2-p) and 2/p, 
\0fAx)\\ P P <c /[.M (/)(*)] (2_P)P/2 - 



{\h\<t>/*} 



e- 2 (f t (x+h)-Mx)f 



dh 



\h\ d +< 



■ dt 



p/2 



m(dx) 



< 



l'[M(f)(x)} p m(dx 



1 (2-p)/2 



e p-2 



< c 



{\h\<fV°} 

\ p P (2 - p)/2 wm 2/2 - 



(ft{x + h) - ft{x)f dtm(dx) 



P/2 



Then the desired result follows. 
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